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Abstract 



^— V ' The linear super-Burnett coefficient gives corrections to the diffusion equation 

^^ . in the form of higher derivatives of the density. Like the diffusion coefficient, it 

^^ I can be expressed in terms of integrals of correlation functions, but involving four 

Qv ' different times. The power law decay of correlations in real gases (with many 

On ■ moving particles) and the random Lorentz gas (with one moving particle and fixed 

^ ', scatterers) are expected to cause the super-Burnett coefficient to diverge in most 

J^' cases. Here we show that the expression for the super-Burnett coefficient of the 

I ■ periodic Lorentz gas converges as a result of exponential decay of correlations and 

Q I a nontrivial cancellation of divergent contributions. 



1 Introduction 

To leading order, the equations of hydrodynamics describe the transport of microscop- 
ically conserved quantities supplemented by linear constitutive relations between the 
relevant thermodynamic forces and fluxes. For example, the conservation of the number 
of particles of a light species moving in an arrangement of fixed scatterers is described 
macroscopically using the continuity equation 

relating the number density n to the current J. Under a hydrodynamic approximation 
where n is "smooth" , hence all of its derivatives are small, the flux J is given as a linear 
function of the force Vn, that is, in an isotropic homogeneous medium, 

J = -DVn (2) 
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leading to the diffusion equation 

DV^n (3) 



dn 2 



dt 

where D is the diffusion coefficient. Similar considerations involving energy and niomen- 
tum conservation lead to the Navier-Stokes equations for a fluid. 

There have been many attempts to go beyond this leading hydrodynamic approxima- 
tion. The constitutive relation can contain an expansion in terms of higher derivatives 
of the forces and/or higher powers of the forces. In the former case one talks of linear 
Burnett [Q] coefficients, and in the latter, of nonlinear Burnett coefficients. The first 
correction in either case (coefficient of V^n or (Vn)^) is often called a Burnett coefficient 
and the second correction (coefficient of V^ra or (Vn)^) a super-Burnett coefficient, and 
so on. We consider the Lorentz gas, in which the light particles do not mutually interact, 
so only linear relations are possible. In addition, it is equally likely for a particle to move 
in one direction or the opposite direction, so the V^n term cancels and we are left with 
the linear super-Burnett coefficient. We will dispense with this cumbersome terminology 
and call it simply the Burnett coefficient. 

Using the Boltzmann equation, hence a limit in which all collisions are uncorrelated, 
the Chapman-Enskog method allows the computation of D as well as other and higher 
order transport coefficients of a low density gas. One of the greatest surprises of nonequi- 
librium kinetic theory was the discovery |^ that density corrections due to recoUisions 
of several particles diverge. This leads to a nonanalytic dependence of transport coeffi- 
cients on density in the case of the random Lorentz gas [|TT| , where noninteracting point 
particles collide with a random arrangement of non-overlapping hard spheres in two or 
three dimensions, but not in an exactly solvable one dimensional model [^ . 

Subsequently the divergences were related to a power law decay of correlations (see [Q 
and references therein), the "long time tails", which are due to density fluctuations in 
infinite random systems. Since transport coefficients can be computed using integrals of 
correlation functions (see below), a slow decay of correlations can lead to divergences. 
As discussed in a review of random Lorentz models [|I|, the diffusion coefficient is typi- 
cally defined, but Burnett coefficients may not be since the latter have more divergent 
expressions in general. The degree of divergence is greater for real gases than for Lorentz 
gases, and decreases with the dimension of the system. Here we show that, as a result of 
the exponential decay of correlations in the periodic Lorentz model, the combination of 
four-time correlation functions used to compute the Burnett coefficient decays sufficiently 
rapidly to ensure convergence. Thus a Burnett coefficient is defined for a random Lorentz 
gas in a dilute approximation and for a periodic Lorentz gas under general conditions. 

2 The Burnett coefficient 

Our notation follows that of Gaspard who gives an account of deterministic diffusion 
including expressions for the Burnett coefficients. He suggests stretched exponential 



bounds on the four time correlations which we find, but does not give the precise form 
or rigorous proof. 

We consider a periodic Lorentz gas in which a point particle moving with unit velocity 
undergoes specular refiections with a periodic array of scatterers. We consider only the 
case of finite horizon, so that time between collisions is bounded. The most usual case is 
that of circular scatterers in a hexagonal lattice, but we demand only that the scatterers 
be convex and sufficiently smooth, specifically C^. We allow for any dimension d > 2, 
but note that a single spherical scatterer per unit cell cannot satisfy the finite horizon 
condition ioT d > 2. Thus we need either nonspherical scatterers or more than one per 
unit cell when d > 2. 

Because the scatterers lie on a periodic lattice, it is possible to represent the dynamics 
in terms of a map (f) acting at the surface of the scatterers. We denote the position and 
velocity direction at the surface of a scatterer in the elementary cell by x. Then the 
lattice translation vector a(x) gives the lattice displacement associated with the free 
fiight following a collision at x; this is a linear combination of the lattice basis vectors 
e*-""-* (a = 1 . . .d) with integer coefficients and is a piecewise constant function of x. T(x) 
is the time for this free fiight. We can compute time averages, which for almost all initial 
conditions are equivalent to the natural measure on the surface of the scatterers () due 
to ergodicity. We define AT(x) = T(x) — (T) so that (AT) = 0. We use Greek indices 
to refer to spatial components such as Oq, and Latin indices to refer to collisions, with 
the shorthand notation aj^ = aa{4>^^) and AT* = AT((/)*x). We have (aj^) = and also 
odd combinations such as J2'ij=-oo{(^a^p'^y) = since time reversibility and phase space 
volume conservation imply that there are equal and opposite contributions from time 
reversed trajectory segments. 

The general periodic Lorentz gas is homogeneous but not isotropic at large scales, 
so the hydrodynamic equation contains transport coefficients that are fully symmetric 
tensors, but do not depend on position: 

OTl 

— = DapVa'^pn + Bap^iVa^pV^Vsn + . . . . (4) 

where B^iS'ys are the Burnett coefficients and repeated Greek indices are summed. A 
general solution can be obtained with a Fourier-Laplace transform of the form n ~ 
exp(st + ik • r) leading to a dispersion relation 

S(k) = -Da(3kak/3 + Baf^^skak/sk^ks + . . . (5) 

The Green-Kubo relations for the transport coefficients in terms of integrals (sums) of 
correlation functions are given in continuous time in, for example, Ref. [Q], and in discrete 
time they are 

\ / i=—oo 
Bafi-yS = — [-^a/375 + 4:C{Dai3D^s + -Da^-D/35 + DasDp^) (7) 



and 



with 



-2{DapE^S + Da-yE/3S + DasE/3^ + Dp^EaS + DpsEa-y + D^sEap)] 



-| oo 

C = -— y {AT'>AT) 



: — CXD 

oo 



Eap = ^^ Y: (AT°a>^^) (9) 



ij=— oo 
oo 



-| oo 



(T) 

-(a>^^)(a>f)-«a^)(a;,a^^) (10) 



ij,fc=— oo 



It is possible to make the lower limits of the sums all zero with some extra symmetrization 
and special treatment of the zero terms. The convergence of the sum for C follows 
in a straightforward manner from the known exponential decay of correlations in this 
system . The convergence of the sum for Eaf3 follows from Theorem 3 in the following 
section. The sum for Fa/B-yS is still less clear since each term diverges; however Theorem 
1 in the following section show that the combination converges, leading to finite Burnett 
coefficients. 

In the case of the hexagonal Lorentz gas, the six-fold rotational symmetry is sufficient 
to restrict the coefficients to their isotropic forms Da/3 = D6a/3 and Ba/3-ys = B^Sa/sS-ys + 
S 0-^^05 + SasSjSy)/^- It seems likely that a straightforward application of our results could 
also show the existence of the higher order Burnett coefficients, for which there exist long 
expressions similar to those given here. In this case it would be interesting to consider 
the radius of convergence of Eq. (|) . 

3 Proofs 

The billiard map (f) is defined on the collision space M of the periodic Lorentz gas. In the 
case of finite horizon, the map (/> and the above functions Oq, and AT are piecewise Holder 
continuous |]^, ^]. This means that M can be partitioned into finitely many domains with 
piecewise smooth boundaries on each of which the function is Holder continuous. 

We always denote by (■) the integration over M with respect to the invariant equi- 
librium measure. In all theorems below, /o,/i,... denote piecewise Holder continuous 
functions on M such that {fk) = for all k. For i > we put fl = fk° 0'- 

Theorem 1 The following series 

oo 

E ifofi^n^m - {fofi'){f2'n') - ifomifi'n') - (/o/^^)(/i7^^) (n) 

always converges absolutely. 



We prove general estimates on multiple correlation functions, from which the above 
theorem will follow. 

Theorem 2 Let ii < ■ ■ ■ < ik O'^d 1 <t <k — 1. Then 

i(/r---/fc'=)-(/r---/r)(i*;r---/fc'=)i<c^fe-Nfc-^irAi^---i^'^ (12) 

where C^ > depends on the functions fi, . . . , fk, and \ < 1 is independent of k and 
/i) • • • ) /fc- 

This allows us to "decouple" multiple correlations provided the gap between the time 
moments it and it+i is long enough. 

Remark. In view of recent advances in the study of correlations for billiard dynamics 
|1^, P] it is very likely that the above estimate can be upgraded to CfcA'**+^~**' for some 



A < 1, which would be an exponential decay of (multiple) correlations. However, the 
estimate (|T^) is easier to prove, and for our purposes it is good enough. 

We postpone the proof of the Theorem ^ and first obtain Theorem |I| and other 
interesting implications. 

Theorem 3 Let ii < ■ ■■ < i^. Consider the k — 1 differences rrit = it+i — it for 
1 < t < k — 1. Let m(i) < m^2) < ■ ■ ■ < '^(fc-i) be the ordered differences mi, . . . , rrik-i- 
Let m = m(j.) where r = [{k + l)/2]. Then 

m■■■fl')\<Ck■\^u-^l?\'^"'' 

where Ck > and A < 1 (here A is the same as in the previous theorem). 

To prove this theorem, we use a simple lemma: 

Lemma 4 Let k > 2 andr = [{k + l)/2]. If we partition k objects into k — r + 1 nonempty 
groups, then at least one group contains exactly one object. 

Now, Theorem ^ is proved by applying Theorem ^ to the k—r largest gaps m(^k-i) , ■ ■ ■ , ''^(r) 
and recalling that (fj) = for every j. 

Note that there is no similar bound on multiple correlations in terms of the next 
difference, m(r+i). Indeed, no matter how large m^r+i) is, we can always arrange the 
time moments ii, . . . ,ik so that ii = i2, is = i^, ■ ■ ., ik-i = "ik, and then the correlation 
(/r ' ' ' Tk) need not be even small, in general it will stay bounded away from zero. 

Corollary 5 For k = 2, we obtain the known j|/ bound on the 2-correlation function: 

mfn<c{f,j2)-\t'-''^''' 

for some C(/i, /2) > and Aq < 1. 



We now prove Theorem |1|. We can assume that ii < 12 < ^3 in (|Il|) and put ffi = 
max{ii,i2 — 'i'l^'is ~ '^2}- Then it easily follows from Theorem ^ and the above corollary 
that 

lif^fi'ftn') - {fon'){f?ft) - {fomifi^m - (/o/r)(/r/^^)i < cm^x^'^' 

This easily implies the absolute convergence of the series (pA|). Theorem [1| is proved. 

Theorem 6 Let Sn = fo + fo + ■ ■ ■ + fo^ ■ Then the k-th moment of Sn can be bounded 
by 

where C(/o, k) > 0. 

This theorem follows from Theorem ^. In the case k = 4 see the proof in 0. The 
argument generalizes to any k > 3 rather directly. In view of the remark made after 
Lemma ^ the above estimate cannot be improved, in general. 

Proof of Theorem |^. We use the partition method developed in and further en- 
hanced in 1^. Note that [|], ^ discuss two-dimensional periodic Lorentz gases, but all 
the results hold in any dimensions as well [^. We recall the essentials of the partition 
method in a more general setting than our Lorentz gas model. 

Let (/) : M — > M be a measurable transformation of a metric space M preserving a 
nonatomic Borel probability measure fi. Let A = {Ai,A2, . . .} be a finite or countable 
measurable partition of M, we denote diam^ = supjjdiamAj}. We put An = 4>^"'A = 
{(f)~"'Ai} for ra > and An,k = An V ■ ■ ■ V Ak ior k > n > 0. A measure of dependence 
between two partitions A = {Ai} and B = {Bj} of the space M is defined to be 

PiA, b) = Y: ha, n B,) - ^l{A,)^,{B,)\ 

hi 
Based on this measure, we put 

/?7v(n) = max /?(A,fc, A+n,JV-i), (13) 

0<fc<A'— n— 1 

for any N >n>{]. Note that l3]\i{n) is called the mixing coefficients of the partition A, 

cf. i- 

For any measurable bounded function f on M and a partition ^ of M we denote 
/^ = (/|v4.) the conditional expectation of / with respect to the cr-algebra generated by 
A, and A^/ = f - f^. We put 

nf{d)= sup IIA^/lb 

Clearly, 'Hf{d) monotonically decreases to zero as c? — *> 0, and the speed of decrease 
characterizes the "smoothness" of /, see P and below. 



Lemma 7 For any functions fi,...,fk on M and any partition A with d = diain^ we 
have 



m ■■■fj:')-{fi'---fh{fi'^--- 11:^)1 <c. 



Y.'HfM + Ph-hi^t+i - it) 



u=i 



(14) 



Proof. First, we replace ft by ft = ft,A- The error can be estimated by 

m ■ ■ ■ fi') - (/r ■ ■ ■ fit)\ < c,Y.\\^Aft\\2 < c.Y.'HfAd) 

t t 

which is a simple calculation using Schwarz' inequality, cf. [Q . Two other correlations in 
(0) can be handled similarly. 

Next, it is rather straightforward that 

K/r---/^'=)-(/r---i*")(it;Y---/fc^)l<nil/*llooxA.-n(^m-^*) 

t 
(for k = 2, this was done in [Q, the general case is similar). Lemma is proved. □ 

We now need a good partition A for our Lorentz gas map (p : M ^ M. Such a 
partition was constructed in P, ^ and further refined in p. Here we state the result: 

Theorem 8 (see ||5|]) Let <p : M -^ M be the billiard map for a periodic Lorentz gas 

with finite horizon. Then for any N > m > 1 there is a partition A = AN,m such that 

(i) diam^ < ci\f^% 

(ii) P^im) < c^N^Xf" 

for some constants Ci, C2 > and Ai, A2 < 1. 

It remains to estimate the term 7if^{d) in (0). Let / be a piecewise Holder continuous 
function on M. Then one can easily obtain |^ that 7if{d) < C{f) ■ c?" for some a > 0. 
Now Theorem ^ follows from Lemma |^ and Theorem ^ 
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